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^ ■ Abstract 

O ' 

CT\ • We calculate the contributions of the axial current to top quark pair production 

in e + e~ annihilation at threshold. The QCD dynamics is taken into account 
by solving the Lippmann-Schwinger equation for the P wave production using 



X 



the QCD potential up to two loops. We demonstrate that the dependence of 



the total and differential cross section on the polarization of the e + and e 
beams allows for an independent extraction of the axial current induced cross 



section. 



Top quark production at an electron-positron collider M has been demonstrated 
to be ideally suited for a precise determination of the top quark mass and for the 
study of its couplings in production and decay. Due to its rapid decay large distance 
nonperturbative QCD effects are irrelevant for the description of the top quark ||, 
and the tt system is well described by perturbative QCD 0. It allows to explore the 
interquark potential at small distances, which is closely related to the strong coupling 
constant. One might eventually even become sensitive to the t-t-Higgs coupling through 
virtual corrections. In order to constrain this multitude of parameters in an optimal 
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way and to reduce inevitable theoretical uncertainties, it is desirable to measure a 
large variety of different observables. Originally the main emphasis had been put on 
the total cross section ||, |J. The excitation curve with its steep rise (the remnant of 
the 15 toponium resonance) is ideally suited for the measurement of the top quark 
mass m t . The correlation between m t and the strength of the potential (a s ) can be 
reduced by comparing data and predictions for the momentum distribution of the top 
quarks ]5], |6|, |7|, [8] , which reflects essentially their Fermi motion in the bound state and 
the smearing of the momentum due to the large decay rate T t , a consequence of the 
uncertainty principle. All these quantities were calculated for the S wave amplitude, 
which is induced by both the electromagnetic current and the vector part of the neutral 



current close to threshold. Expanding in the limit of small velocities (3 = y 1 — 4mj/s 
(\/s being the total centre of mass energy), the next term is due to S — P wave 
interference. The subleading P wave amplitude originates from the production through 
the axial vector current. The interference term is responsible for the anisotropic angular 
dependence, specifically the term linear in cos 6, and the resulting forward-backward 
asymmetry ||. Similarly, an angular dependent polarization of top quarks is induced 
by the S — P wave interference which adds to the dominant polarization parallel to the 
e + e~ beams []10| . Rescattering corrections [TT], , although important for the detailed 
quantitative analysis, do not alter this qualitative picture. 

Clearly, the next step in this sequence of improvements are corrections of order f3 2 
which, for interacting quarks close to threshold, translate into corrections of order a 2 
and (3a s . For the vector current this has been recently persued by different groups, 
which have demonstrated the importance of these next-to-next-to leading order cor- 
rections |I3, [14], [Tj|. However, in the same order (3 2 (or a 2 ) also axial vector induced 
contributions must be incorporated. They affect both the excitation curve and the 
momentum distribution. Close to threshold these axial contributions are suppressed 
relative to the dominant 5* waves by two powers of (3 whence a treatment of the leading 
terms is sufficient for the present purpose. These axial contributions are mediated by 
the virtual Z boson only. Therefore their dependence on the beam polarization differs 
from the one of the vector current induced rate. This, in turn, allows for the separation 
of the two independent contributions to the total and differential cross section. With 
the axial contribution representing an independent observable, this separation is pos- 
sible independent of potential uncertainties in the NNLO calculation of the dominant 
piece. However, in view of the f3 2 suppression of the axial rate and the relatively small 
couplings of the neutral current, large luminosities and a high degree of polarization 
are required to make a clean extraction of the axial part possible. These features are 
unique for linear colliders, as proposed e.g. in [[!], [T6| . However, even without this pos- 



sibility, it will be important to control the impact of this contribution on the extraction 
of the top quark mass and the interquark potential. Let us also stress that the axial 
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rate, although closely related to the S — P wave interference piece, is an independent 
observable. Rescattering corrections, which are present in the angular distribution and 
in the top polarization, are calculated to 0(a s ) [11], [12] but shown to be unimpor- 



tant as long as the total cross section is concerned |T7|, [18|]. In addition, rescattering 
corrections do not affect the separation between axial and vector contributions. 

P wave threshold production of massive quarks in 77 collisions has been analysed for 
the case of a pure Coulomb potential in Ref. |19| and much of the general considerations 
can be taken over to the present case.Q This refers in particular to the treatment of the 
linearly divergent integrals over the momentum distribution and the order of magnitude 
estimates. However, for definite predictions the QCD potential with its logarithmically 
varying coupling strength has to be employed. The relative size of the electromagnetic 
and weak couplings is important for the phenomenological analysis, as well as the 
dependence on the beam polarization. 

The momentum distribution of the top quark including the influence of beam po- 
larization can be written in the form 



da 
dp 



ml 



P+P-) [(ax + xa 2 ) (l - ~) Vss(p, E) 



where the correction factors from hard gluon exchange, (1 — 16a s /37r) and (1 — 8a s /37r), 
are taken from |]2"2"| , |2"B| | . P + and P denote the polarization of the positron and electron 
beams, respectively, and x is defined as 



X 



P+-P- 
1 - P+P_ 



The coefficients a,- read 



a i = {QeQt + v e v t df + (a e v t dY 
a 5 = (a t df (yl + afj , 



a 2 = 2a e v t d (q e q t + v e v t d) 
a 6 = 2v e a e (a t d) 2 , 



with 



d 



16sin 2 ^cos 2 ^ s -Ml 
and the electromagnetic and weak charges 



(3) 
(4) 



qt 



■--I, 

2/3, 



v e 



-1 +4 sin 2 6 W , 
1 - 8/3 sin 2 W 



(h 



1 . 



(5) 



1 See also Rcfs. p(j and [EE for related discussions of P wave production of quarks and squarks in 



collisions near threshold. 



3 



The dynamics of the strong interaction is contained in the functions 



V S -s(p, E) = p 2 \G(p, E)\ 2 and V P _ P {p, E) = p 



— Hp, e) 

m t 



(6) 



where E = y/s — 2m t is the energy relative to the nominal threshold. The S and P 
wave Green functions G(p, E) and F(p, E) fulfill the Lippmann-Schwinger equations 

G(p, E) = G (p, E) + G (p, E)J^-V{ \ p- q\ ) G(q, E) , (7) 
F(p,E) = G o (p,E)+G o (p,E)J-0- 3 ^-V(\p-q\)F(q,E) (8) 

where p = \p\ is the momentum of the top quark in tt rest frame, Gq(p,E) = 
( E — p 2 /m t + iT t ) 1 is the free Green function, and T t denotes the top quark width. 
For the QCD potential in momentum space, V, we adopt the two loop result |24| with 
the long distance regularization as described in [ 25| . Eqs. (0, are then solved numer- 



ically as described in [EL For all the results discussed below we use the parameters 
m t = 175 GeV, r t = 1.43 GeV and a s (M 2 z ) = 0.118. 

For large momenta both G(p, E) and F(p, E) approach the free Green function Go- 
It is thus evident that the integral over the momentum distribution diverges linearly 
for the P wave. This is, however, an artefact of the nonrelativistic approximation. 
The problem could be cured, for example, by introducing in this region the relativistic 
(free) Green function and phase space and by treating the interaction as a (small) 
perturbation. However, in practice, a cutoff will be provided by the experimental 
analysis. The invariant mass of the W plus b jet in events with large p (Wb) and small 
E = \fs — 2m t will necessarily be strongly shifted away from m t towards smaller values. 
Such events will either not be included in the ti sample or, in any case, will require 
special treatment. Hence, wherever total cross sections are presented, a cutoff p max of 
order m t will be introduced which is easily included also in the experimental analysis. 

The relative magnitude of the P wave result is best visualized by considering the ba- 
sic elements T>s-s and T> P _ P which enter Eq. ([]]). In Fig. |l] we show these distributions 
for three energies, E = —3, and 3 GeV. These energies roughly correspond to the 
location of the 15 peak, the nominal threshold and the onset of the continuum. For the 
S wave (Fig. |l]a) we observe a fairly wide distribution at E = — 3 GeV, a consequence 
of the momentum spread of the constituents in the IS bound state. With increasing 
energy the interaction becomes less important, the width of the distribution decreases 



and approaches the free result T t ym t /E. For the P wave (Fig. |T]b) the contribution is 
tiny at E = —3 GeV and develops a peak only gradually with increasing energy. The 
ratio T>s-s/'Dp-p is shown in Fig. |l]c. For energies well above threshold its behaviour 
is essentially given by the factor p 2 /m 2 , since both F and G are approximated by the 
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free Green function Go- However, for energies relatively close to threshold the strong 
interaction distorts the free wave functions which leads to a deviation from the pure 
p 2 /m 2 behaviour. The ratio of the integrated S and P wave distributions as functions 
of E is shown in Fig. 0a. The different curves give the results for different values of 
the momentum cutoff p max which is applied both in numerator and denominator. For 
a realistic analysis p max = m t /3 or m t /2 should be used at most. For free and stable 
quarks the ratio is given by p 2 (E)/m 2 ~ E/m t . Fig. |2]a shows that close to threshold 
the momentum spread from the QCD bound state dynamics leads to a significant mod- 
ification of the E/m t behaviour and increases the P wave contribution. The minimum 
of the ratio Jo' max dpV P _ P / Jo >ma3t dp V S -s is observed roughly at the location of the 
remnant of the IS" peak of the R t i ratio (Fig. |||b) 



'max 



Rtt=—r- " r = 2 / dpVss- 9 

a(e + e~ — > 7* — ► /v/i ) Timf Jo 

With these ingredients it is now straightforward to predict the differential distribution 
for the three characteristic polarizations P = —1, 0, +1 and P + = 0. The cross 
sections are drastically different for the three choices, see Fig. |3], reflecting the large 
left-right asymmetry = a 2 /ai ~ 0.4 of the S wave contribution 0, |27J. Including 



the small P wave term (dotted curves) leads to marginal changes only, which are barely 
visible in Fig. § even for the highest energies. The relative size of the axial contribution 
is better visible in Fig. f| where the ratio between the axial and the vector contribution 
is plotted as a function of the momentum p. The shapes and the magnitude are 
fairly similar for the different energies. This is a consequence of the fact, that the 
ratio Up_ P (p, E) /V S -s(p, E) is relatively insensitive to the energy. In fact, in the 
absence of interaction this ratio is just given by p 2 /m 2 , independent of E. In contrast, 
the location of the maximun of the distribution itself varies with E, and this is mainly 
responsible for the increase of the integrated P wave cross section. The integrated cross 
section with and without the P wave contribution is shown in Fig. [5]a, where for the 
cutoff p m ax = w*t/2 is adopted. The ratio between axial and vector contributions, both 
integrated up to m t /2 is shown in Fig. |5|b. The shape of these curves reflects the shape 
of the ratio / dpVp^p / J dpVs-s displayed already in Fig. |2|a. The normalization 
depends on the polarization. This demonstrates that experiments with polarized beams 
are able to extract <j££ separately, provided that a statistical and systematic precision 
at the percent level can be reached. In any case, if a theoretical prediction of shape 
and normalization of da /dp at a precision of one or two percent is needed the P wave 
contribution to the cross section should be included. 
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Figure 1: Results for the basic elements in Eq. ([!]): a) Vs-s(p, E), b) T>p-p(p, E) and 
c) the ratio Vp_p /T>s-s for the three energies E = — 3 GeV (continuous curves), E = 
(dashed lines) and E = 3 GeV (dotted) as a function of the top quark momentum p. 
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Figure 2: 

a) Ratio of the integrated distributions: Jo max dpV P _ P (p, E) / Jo max dpT>ss(p, E) as a 
function of the energy E = yfs — 2m t for four different values of the cutoff: continuous, 
dashed, dotted, and dash-dotted lines correspond to p max = mt • |, \, |> 1 , respec- 
tively, b) The normalized total cross section R t i as defined in Eq. @) as a function of 
E. 
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Figure 3: Differential cross section da(e + e~ — > ti)/dp as defined in Eq. ([!]) as a 
function of p for six different energies, E = —3, 0, 3, 5, 10, 20 GeV, as indicated in 
the plots a). . . f). The continuous, dashed and dash-dotted lines show the pure S 
wave result for the three different choices of the e~ polarization P_ = — 1, and 1, 
respectively. (P + = 0.) The dotted lines show the full result including the P wave 
contributions. 





Figure 4: Relative size of the axial contribution da AA /dp compared to the vector 
contribution da vv /dp to the differential cross section as a function of p for six different 
values of the energy E [plots a). . . f)]. The continuous, dashed and dash-dotted lines 
correspond to P_ = — f, and 1, respectively. (P+ = 0.) 
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continuous, dashed, dash-dotted: P = -1, 0, +1 (P + = 0) 
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Figure 5: a) The total cross section a(e + e~ — > tt) as a function of E for three 
different choices of the e~ polarization: the continuous, dashed and dash-dotted lines 
correspond to P_ = — 1, and 1, respectively, where only S wave production is taken 
into account. The dotted lines show the corresponding total cross sections including 
the P wave contributions, b) Ratio of the P to the S wave contribution a^/a^ for 
the three different e" polarizations. 
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